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In  linear  regression,  gross  rainfall  alone  accounted  for  85  percent  to  as  high  as  99 
percent  of  throughfall  variation  at  individual  gage  points.  Stemflow  was  also  primarily 
dependent  on  gross  rainfall.  Results  further  demonstrate  adjustment  of  throughfall 
for  mean  canopy  density,  adjustment  of  stemflow  for  tree  size  (d.b.h.),  and  combination 
of  these  relationships  to  estimate  net  rainfall  for  different  stand  densities. 
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Variabilities  of  throughfall  and  stemflow 
have  been  variously  attributed  to  random  error 
(Stout  and  McMahon  1961),  tree  parameters, 
or  characteristics  of  the  forest  canopy.  In 
their  extensive  review  of  interception  studies 
in  mature  hardwoods,  Helvey  and  Patric  (1965a) 
concluded  there  is  no  consistent  evidence  that 
interception  losses  are  greatly  affected  by  a 
variety  of  canopy  densities.  They  go  on  to 
say,  however,  that  failure  to  show  variations 
in  throughfall  has  probably  been  due  more  to 
failure  of  sampling  methods  to  adequately  meas- 
ure variation  than  to  lack  of  variation  itself. 
They  further  maintain  that  throughfall  under 
mature  hardwoods  must  be  inversely  propor- 
tional to  canopy  density.  The  same  logic  must 
also  apply  to  conifers. 

Stout  and  McMahon  (1961)  conclude  that 
amount  of  throughfall  under  a  specific  tree 
may  vary  with  position  or  direction  of  the 

1  Hydrologist ,  Rocky  Mountain  Forest  and  Range 
Experiment  Station,  located  at  Rapid  City,  South  Dakota, 
in  cooperation  with  South  Dakota  School  of  Mines  and 
Technology;  Station's  central  headquarters  maintained 
at  Fort  Collins,  in  cooperation  with  Colorado  State  Uni- 
versity. 


sampling  point  from  the  trunk.  Their  findings 
agree  in  general  with  earlier  works  reported  by 
Beall  (1934),  Ovington  (1954),  and Geiger  (1957). 
Such  variations  have  been  attributed  to  density 
of  foliage  immediately  above  the  gage  (Horton 
1919).  Most  such  conclusions,  though  logical, 
are  not  well  supported,  and,  as  pointed  out  by 
Helvey  and  Patric  (1965a),  do  not  adequately 
define  relationships  of  an  interlaced  forest  cano- 
py, especially  where  there  is  more  than  one 
canopy  level. 

Stemflow  is,  as  a  rule,  more  variable  than 
throughfall.  Helvey  and  Patric  (1965b)  attribute 
this  greater  variability  to  the  endless  variety 
of  branch  arrangements  and  bark  roughness 
between  trees.  In  some  studies  there  have 
been  no  strong  correlations  of  several  tree 
dimensions  with  stemflow  (Black  1957).  In  a 
number  of  other  studies  summarized  by  Helvey 
and  Patric  (1965a),  stemflow  has  been  definitely 
associated  with  such  factors  as  bark  roughness, 
trunk  diameter,  and  height  of  crown  above  the 
general  canopy  level.  Lawson  (1967)  concluded 
that  addition  of  a  tree  size  variable  greatly 
improved  stemflow  predictions  for  both  pines 
and  hardwoods. 
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Opportunity  for  further  analysis  of  the  in- 
fluence of  canopy  and  tree  variables  on  through- 
fall  and  stemflow  was  provided  by  measurements 
taken  in  connection  with  a  study  of  soil  mois- 
ture trends  after  clearcutting  and  thinning  in  a 
second-growth  ponderosa  pine  stand  in  the 
Black  Hills  (Orr  1968).  The  accumulating 
evidence  that  interception  reduces  transpiration 
to  some  degree  (Rutter  1968)  makes  it  even 
more  important  that  the  throughfall  and  stem- 
flow  processes  be  understood  and  accounted 
for  in  hydrologic  analyses. 


Study  Area 

Throughfall  was  measured  on  three  30-  by 
30-foot  subplots  each  in  a  120-  by  150-foot 
thinned  plot  and  in  an  adjacent  unthinned  plot 
in  second-growth  ponderosa  pine  about  70 years 
old.  The  thinned  stand  had  435  trees  per 
acre  averaging  36  feet  tall  and  5.8  inches  diam- 
eter, breast  high  (d.b.h.).  The  adjacent  un- 
thinned stand  contained  2,885  trees  per  acre 
averaging  29  feet  tall  and  3.5  inches  d.b.h. 


Stemflow  was  measured  on  all  10  trees  on 
one  30-  by  30-foot  thinned  subplot,  and  on  a 
random  sample  of  10  out  of  54  trees  on  one 
30-  by  30-foot  unthinned  subplot. 

The  plots  were  on  the  lower  and  more 
gently  sloping  (9  percent)  portion  of  a  steep 
north-  to  northeast-facing  slope.  Soils  appeared 
to  have  developed  in  place  from  limestone 
parent  material.  Site  index  is  between  70  and 
75,  near  maximum  for  ponderosa  pine  in  the 
Black  Hills.  Precipitation  averaged  21.3  inches 
per  year  over  a  5-year  period.  Precipitation 
usually  is  minimum  during  the  winter,  builds 
to  a  maximum  in  June,  and  declines  gradually 
through  the  summer  and  fall. 


Methods 

Throughfall  was  measured  with  two  standard 
8-inch  cans  rotated  about  a  system  of  10  ran- 
dom points  on  each  of  the  three  subplots  in 
each  of  the  two  treatments.  Hence,  30  points 
were  sampled  in  each  treatment.  Measure- 
ments were  taken  yearlong  at  minimum  1-week 


Two  standard  8-inch  cans  were  rotated 
about  a  system  of  10  random  points  on 
each  of  three  subplots  in  each  of  two 
treatments  to  measure  throughfall . 


intervals  for  3  years.  Occasional  in-between 
measurements  were  taken  in  case  of  large 
storms.  Gages  were  weighed  and  contents 
converted  to  inches  depth.  The  two  gages 
on  each  subplot  were  moved  to  another  pair 
of  points  after  any  measurement  period  in 
which  gross  precipitation  equaled  or  exceeded 
0.05  inch.  Because  gages  were  located  at  only 
two  of  10  sampling  points  per  subplot  at  any 
given  time,  there  are  five  different  sets  of 
precipitation  data,  each  representing  a  different 
set  of  six  gage  points. 

Stemflow  was  collected  in  copper  collars 
attached  to  the  trees  and  sealed  with  plastic 
cement.  Collar  openings  were  about  one-half 
inch.  The  water  was  piped  to  6-gallon  con- 
tainers and  weighed  for  each  tree.  Stemflow 
was  measured  for  2  years. 

Gross  precipitation  was  measured  with  both 
a  recording  and  standard  gage  in  a  nearby 
forest  opening  with  minimum  45°  clearance, 
and  in  one  standard  gage  near  the  center  of  a 
120-  by  150-foot  clearcut  plot  adjacent  to  the 
thinned  and  unthinned  plots.  Oil  was  used  in 
all  gages  to  minimize  evaporation. 

Canopy  was  photographed  from  ground  level 
at  each  of  the  30  throughfall  gage  points  in 
each  of  the  two  treatments,  with  the  hemi- 
spherical camera  described  by  Brown  (1962). 
Percentage  canopy  was  estimated  in  each  of 
80  equal  hemisphere  segments  from  zenith  to 
90°.  Values  were  averaged  for  an  assortment 
of  zenith  and  azimuth  segments.  Estimates 
were  made  by  two  different  individuals  and 
repeated  by  one  individual  after  several  months 
had  elapsed.  There  was  practically  no  difference 
in  estimates  between  individuals  or  between 
repeat  estimates  by  the  same  individual. 

Tree  measurements,  including  d.b.h.,  height 
to  live  crown,  total  height,  and  crown  diameter 
were  made  by  standard  mensurational  tech- 
niques. 

Only  the  summer  rainfall  data  (May-Septem- 
ber, 1958-60  inclusive)  were  analyzed  in  detail. 
Winter  and  spring  data  were  complicated  by 
occasional  mixed  rain  and  snow,  and  by  re- 
cording difficulties. 

The  principal  throughfall  analyses  involved 
stepwise  regression  of  throughfall  on  rainfall 
variables  for  each  gage  point,  followed  by 
regression  of  mean  throughfall  at  individual 
points  on  canopy  density.  The  precipitation 
variables  tested  in  the  initial  analyses  were 
gross  rainfall,  number  of  storms  per  observation 
interval,  and  an  expression  of  rainfall  intensity. 
In  the  final  analysis,  all  30  sampling  points 
in  each  treatment  and  the  two  treatments  were 
combined. 


Stemflow  analyses  involved  regression  of 
stemflow  on  rainfall  variables  for  individual 
sample  trees,  introduction  of  tree  variables, 
and  final  combination  of  all  sample  trees  in 
the  two  treatments. 

Assumptions  of  linearity  were  checked.  Re- 
sults did  not  justify  an  attempt  to  fit  curves 
in  either  the  throughfall  or  stemflow  analyses. 


Results 

Throughfall 

Rain  throughfall  was  more  variable  in  the 
unthinned  than  in  the  thinned  stand.  Mean 
canopy  density  in  full  azimuth,  52°  zenith  pro- 
jection, ranged  from  28  to  48  percent  and  aver- 
aged 42  percent  in  the  thinned  stand,  and 
ranged  from  52  to  75  percent  and  averaged 
66  percent  in  the  unthinned.  Part  of  the 
greater  throughfall  variability  in  the  unthinned 
stand  was  obviously  due  to  greater  magnitude 
of  sampled  drip  concentration  there.  At  one 
point,  for  example,  where  measured  canopy 
density  was  70  percent,  the  recorded  throughfall 
was  two  to  three  times  greater  than  gross 
rainfall,  but  only  for  gross  rainfall  amounts 
exceeding  about  2  inches.  Similar  results  have 
been  reported  by  other  investigators  (Ovington 
1954,  Wicht  1941). 

In  all  of  the  30  individual  sampling-point 
regressions  in  each  of  the  two  treatments, 
gross  rainfall  alone  accounted  for  92.9  to  99.7 
percent  of  throughfall  variation  in  the  thinned 
stand  and  from  85.5  to  99.7  percent  in  the 
unthinned.  Number  of  observations  per  gage 
point  ranged  from  8  to  12.  The  rainfall  in- 
tensity variable  was  also  highly  correlated  with 
throughfall  at  24  to  30  points  in  both  the 
thinned  and  unthinned  stands,  and  the  number 
of  storms  per  observation  interval  (3-hour 
separation)  was  significantly  correlated  with 
throughfall  at  more  than  half  the  sampling  points 
in  both  treatments.  However,  neither  of  these 
variables  accounted  for  significant  throughfall 
variation  after  gross  rainfall  in  either  treat- 
ment. Canopy  density  did  not  account  for  a 
consistently  significant  proportion  of  through- 
fall  variation  for  different  storm  observation 
groups  in  the  separate  treatments.  The  range 
of  canopy  density  was  apparently  too  narrow. 
When  the  two  treatments  were  combined,  how- 
ever, the  range  was  broadened  sufficiently  for 
valid  expression  of  the  effect  of  canopy  on 
throughfall. 
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The  final  regressions  of  throughfall  (Y]  )  on 
gross  rainfall  (X  ] )  alone  were: 


Thinned 

Y,  =  -0.004  +  0.888X, 
Unthinned 


Y,  =  -0.054  +  0.813X, 


(R    =  0.96)  [1] 
(R2=  0.82)[2] 


Combined  regression  for  thinned  and  un- 
thinned, incorporating  the  52°  zenith  projection 
full  azimuth  canopy  density  variable  [X2]  ex- 
pressed in  percent  yielded  the  following 
equation: 


Y2=  0.167 


0.851X  !  -  0.0037X2  (R  =  0.89)  [3] 


Throughfall  estimates  from  separate  equation 
[1]  and  [2]  are  compared  with  estimates  for 
both  treatments  obtained  with  equation  [3] 
(table  1).  Equation  [3]  estimates  mean  through- 
fall  from  the  mean  of  gross  rainfall  with  adjust- 
ment for  canopy  density. 


Mean  canopy  densities  for  other  zenith  pro- 
jections and  in  the  azimuth  range  of  prevailing 
winds  were  computed,  but  those  tested  were 
interrelated  to  so  high  a  degree  that  no  single 
one  was  significantly  better  than  another.  This 
does  not  mean  that  closer  relationships  or  better 
expressions  of  canopy  density  do  not  exist. 
A  more  efficient  sampling  scheme  specifically 
designed  for  the  purpose  will  be  necessary  to 
establish  such  relationships.  Nevertheless,  the 
present  study  clearly  demonstrates  use  of  canopy 
density  parameters  to  quantitatively  account 
for  a  portion  of  throughfall  variation. 


Stemflow 

All  completely  recorded  amounts  of  gross 
rainfall  larger  than  0.2  inch  yielded  measurable 
stemflow  from  at  least  one  gaged  tree  in  one 
or  the  other  of  the  two  treatments.  The  final 
analysis  involved  21  summer  rainfall  events 
larger  than  0.2  inch  in  both  of  the  treatments 
over  the  2  years  of  measurement. 


Table  1 . --Reg ress ion  estimates  of  mean  throughfall   using  equations  1,  2,  and  3 


I  tern 


Observation  group1 


x2 


Mean  observed  gross  rainfall  (inches) 

Observed  mean  canopy  density  (percent) 
52    zenith  projection,  full  azimuth 

Thinned  stand 

Unthinned  stand 

Mean  observed  throughfall  (inches) 
Thinned  stand 
Unthinned  stand 


Yj     Estimated  mean  throughfall  (inches) 
Thinned  stand  (Eq.  1) 
Unthinned  stand  (Eq.  2) 

Y2    Estimated  mean  throughfall  (inches) 
from  mean  gross  rainfall  with  adjust- 
ment for  canopy  density  (Eq.  3) 

Thinned  stand 

Unthinned  stand 


0.680 


41  .8 
67.3 


59 

.47 


60 
50 


.59 
.50 


0.430 


41 .8 
64.2 


41 
31 


38 
30 


38 
30 


0.751 


41  .0 
66.3 


.66 
.54 


,66 
56 


,65 
56 


0.866 


40.3 
69.0 


77 
58 


.77 
.65 


•  75 
.65 


0.743 


43.3 
61 .8 


63 

,64 


66 
55 


.64 
.57 


Each  group  represents  a  set  of  six  gages,  each  set  having  a  different  group  of  rainfall  events 
due  to  gage  rotation. 
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Here,  as  in  the  case  of  throughfall,  gross 
rainfall  was  the  only  significant  precipitation 
variable  in  combined  regressions.  The  com- 
bination equations  for  regression  of  pounds  of 
stemflow  (Y)  on  gross  rainfall  (X^  for  all  10 
trees  in  each  of  the  two  treatments  are: 

Thinned 

Y  =  -1.94+13.21X1  (R2=0.57)  [4] 
Unthinned 

Y  =  -2.43  +  13.09X1  (R2=0.31)  [5] 

X,  >  0.2  inch 

The  addition  of  tree  variables  greatly  im- 
proved on  these  regressions.  All  three  variables 
tested— d.b.h.,  tree  height,  and  crown  volume 
(taken  as  volume  of  cylinder  whose  diameter  = 
mean  crown  diameter)— were  significantly  cor- 
related with  stemflow,  but  height  and  crown 
volume  were  not  significant  after  d.b.h.  because 
of  high  intercorrelation.  R2  values  increased 
to  0.66  and  0.50  for  thinned  and  unthinned 
stands,  respectively,  after  including  d.b.h. 
D.b.h.  appears  to  be  more  important  in  the 
unthinned  stand,  considering  the  greater  in- 
crease in  R2 . 

The  regression  equation  after  combining 
the  two  treatments  and  incorporating  d.b.h. 
(X2)is: 

Y  =  -21.18+ 13. 10X,  +4.49X2  (R2  =0.55)  [6] 
X,  >  0.2  inch 

Estimated  equivalent  areal  depth  would  be 
equal  to  estimated  pounds  for  the  average  size 
tree  on  a  given  sample  area,  times  the  number 
of  trees,  converted  to  volume  and  divided  by 
land  area.  Size  of  trees  involved  in  this  study 
ranged  from  2.9  to  6.6  inches  d.b.h. 


Throughfall  and  Stemflow 
Combined  (Net  Rainfall) 

Addition  of  throughfall  and  stemflow  yields 
net  rainfall,  here  defined  as  depth  of  rainfall 
reaching  the  surface  of  the  forest  floor.  In 
the  study  area  there  was  virtually  no  under- 
story  vegetation— just  a  bare  mat  of  pine  needles 
at  the  forest  floor  surface.    Combination  of 


equations  [3]  and  [6]  yielded  the  following 
equation  for  net  rainfall: 

Y  =  n,  (0.167 +  0.851X!  -  0.0037X2)  [7] 

+  n2N(-0.0000936  +  0.00005783X3 
+  0.00001986X4) 

where 

Y  =  Net  rainfall  (inches) 

n-\    =  Number  of  recorded  precipitation  events 

>  0.05  inch 

Xj  =  Average  depth  (inches)  of  all  precipitation 

events  >  0.05  inch 
X2  =  Average  canopy  density  (percent) 
n  2  =  Number  of  recorded  precipitation  events 

>  0.20  inch 

X3  =  Average  depth  (inches)  of  rainfall  events 

>  0.20  inch 

X4  =  Average  d.b.h.  (inches)  of  all  trees  on 

sampled  area 
N    =  Number  of  trees  per  acre 


The  first  part  of  the  equation  provides  an 
estimate  of  throughfall  adjusted  for  canopy 
density,  and  the  second  part  an  estimate  of 
pounds  of  stemflow  converted  to  inches  depth 
according  to  tree  d.b.h.  and  number  of  trees 
per  acre.  This  relationship  is  based  on  and 
hence  is  strictly  applicable  to  a  rather  narrow 
situation— one  forest  type,  one  age  class,  and 
one  climatic  regime— in  which  only  precipitation 
amounts  equal  to  or  larger  than  0.05  inch  were 
considered. 

However,  considering  the  surprising  simi- 
larity of  results  from  studies  of  different  species 
at  widely  separated  locations,  it  seems  likely 
that  the  equation  may  yield  realistic  estimates 
of  the  relative  magnitude  of  net  rainfall  (or 
interception  loss)  over  a  much  larger  area. 
For  such  use  the  exact  method  of  determining 
canopy  density  probably  is  not  critical  so  long 
as  it  is  consistent. 


Summary  and  Discussion 

This  study  demonstrates  a  possible  method 
of  estimating  net  rainfall  in  relation  to  stand 
density  of  second-growth  ponderosa  pine  in 
the  Black  Hills.  As  in  virtually  all  reported 
studies  in  both  conifer  and  deciduous  forest 
types,  gross  precipitation  depth  is  the  primary 
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Canopy  was  photographed  from  ground  level  at 
each  of  30  throughfall  gage  points  in  each  of 
two  treatments ,  with  a  hemispherical  camera 
(Brown  1962) : 


This  thinned  stand  has  435  trees  per  acre 
(N)  ,  averaging  5.8  inches  d.b.h.    (x^)  • 
Canopy  density  (xz)  -*s  about  41  percent 
(in  52°  zenith  projection) .     Consider  a 
single  storm  of  0.75  inch   (n\  =  1,  n^  =  1, 
Xi  =  0.75,   X3  -  0.75).     Applying  equation 
[7],  throughfall  is  calculated  as  0.65 
inch  and  stemflow  0.03  inch,  for  a  total 
net  rainfall  of  0.68  inch,  or  91  percent 
of  the  amount  of  rainfall  reaching  the 
ground  in  the  open. 


This  unthinned  stand  has  2,885  trees  per 
acre   (N) ,  averaging  3.5  inches  d.b.h.  (Xk) 
Canopy  density   (xz)  i-s  about  66  percent 
(in  52°  zenith  projection)  .     Consider  a 
single  storm  of  0.75  inch   (n\  =  1,  n2  =  1, 
Xi  =  0.75,   X3  =  0.75).     Applying  equation 
[7],  throughfall  is  calculated  as  0.56 
inch  and  stemflow  0.06  inch,  for  a  total 
net  rainfall  of  0.62  inch,  or  83  percent 
of  rainfall  reaching  the  ground  in  the 
open . 
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controlling  variable  in  both  throughfall  and 
stemflow.  This  primary  control  is  clearly  ob- 
vious in  studies  that  have  involved  regression 
analyses.  However,  a  combination  of  results 
from  adjacent  thinned  and  unthinned  plots 
accounted  for  an  additional  small  but  neverthe- 
less significant  proportion  of  both  throughfall 
and  stemflow  variances .  This  provided  the  basis 
for  adjusting  throughfall  for  percent  canopy 
density  and  adjusting  stemflow  for  tree  d.b.h. 
The  combination  of  these  two  relationships 
yields  an  equation  for  net  rainfall. 

Canopy  density  is  the  most  obvious  factor 
that  it  would  be  expected  might  influence 
throughfall.  However,  other  researchers,  as 
pointed  out  earlier,  have  concluded  that  there 
is  no  consistent  evidence  that  interception  losses 
are  greatly  affected  by  a  variety  of  canopy 
densities.  The  overpowering  influence  of  gross 
precipitation  in  the  ordinary  regression  approach 
is  very  likely  one  of  the  main  reasons  for  this 
lack  of  consistency.  In  the  present  study, 
using  stepwise  regression,  canopy  density  also 
was  not  significant  until  results  from  adjacent 
thinned  and  unthinned  plots  were  combined. 
The  combination  resulted  in  a  broad  enough 
range  of  canopy  density  to  define  a  significant 
relationship. 

A  variety  of  canopy  measurements  involving 
average  percent  density  in  different  zenith  pro- 
jections and/or  azimuth  segments  were  tested. 
Densities  were  estimated  from  vertical  photos 
taken  from  the  ground  upward  at  each  of  the 
throughfall  points,  60  in  all.  Because  of  inter- 
correlation  and  small  residual  variance  after 
gross  rainfall,  no  one  canopy  variable  tested 
significantly  better  than  another.  On  a  rational 
basis  the  average  density  in  52°  zenith  pro- 
jection, full  azimuth,  was  used  in  final  analysis. 

Similar  tendencies  were  evident  in  stemflow 
regression  analyses.  A  variety  of  tree  size 
variables  were  tested,  including  d.b.h.,  total 
height,  and  crown  volume.  Each  correlated 
with  stemflow  independently,  but  in  stepwise 
regression  neither  height  nor  crown  volume 
were  significant  after  d.b.h.  because  of  inter- 
correlation. 

The  foregoing  indicates  statistically  defin- 
able relationships  of  canopy  density  with 
throughfall  and  tree  size  with  stemflow.  Where 
or  when  more  detailed  information  is  needed 
for  solution  of  specific  hydrologic  problems, 
other  sampling  design  and  analysis  techniques 
may  yield  better  defined  relationships.  In  the 
meantime,  entry  of  appropriate  stand  measure- 
ments in  the  combined  equation  [7]  will  yield 
a  realistic  idea  of  the  magnitude  of  net  rainfall 
(or  interception)  in  relation  to  stand  density 
of  second-growth  ponaerosa  pine. 
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